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Fe-P  alloys  with  high  phosphorous  content  have  been  targeted  as  promising  anode  materials  because  of 
their  high  theoretical  capacity.  However,  the  synthesis  and  cycling  performance  remain  great  challenges. 
Hereby  FePy  (3  <  y  <  4)  nanoparticles  are  facilely  synthesized  through  a  dry  mechanochemical  method 
by  reacting  iron  and  red  phosphorus  powders  in  an  inert  atmosphere.  The  morphology  and  crystal 
structure  of  this  material  are  characterized  by  SEM  and  XRD,  respectively,  while  the  electrochemical 
performance  is  evaluated  by  a  number  of  different  techniques.  The  1st  and  2nd  discharge  capacity  of  FePy 
reaches  1984  mAh  g-1  and  1486  mAh  g-1,  respectively,  and  after  10  cycles  at  0.03  mA  cnrr2  (20  mA  g-1, 
0.03C),  the  capacity  remains  1089  mAh  g-1  with  a  coulombic  efficiency  of  97%,  much  higher  than  the 
reported  results  to  date.  The  cyclability  of  this  material  becomes  fairly  better  at  a  higher  current  density, 
but  the  specific  capacity  is  lower  compared  to  that  of  the  smaller  current  density.  By  adding  fluoro- 
ethylene  carbonate  (FEC)  to  the  electrolyte,  the  cycling  performance  of  this  material  was  improved.  The 
EIS  analysis  has  also  been  performed  in  order  to  better  understand  the  capacity  fade  mechanism  of  FePy. 
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1.  Introduction 

With  the  demand  of  portable  electronic  devices  and  electric 
vehicles  growing,  lithium-ion  batteries  with  a  higher  energy  den¬ 
sity  are  attracting  more  and  more  attentions  from  both  the 
fundamental  studies  and  the  real-world  applications.  In  addition  to 
the  electrolyte  and  separator,  the  research  community  has  mainly 
focused  on  the  novel  electrode  materials  with  higher  capacity  and 
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energy  density  to  satisfy  the  demand.  Among  these  electrode  ma¬ 
terials,  alloys  and  intermetallic  compounds  are  promising  candi¬ 
dates  for  energy  storage  owing  to  their  remarkable  advantages 
[1-5  .  Regarding  a  number  of  different  methods  to  synthesize 
these  alloys  and  intermetallic  compounds,  the  mechanochemical 
method  is  the  first  choice  and  is  always  being  utilized  to  prepare 
nanoscale  alloys  because  it's  facile,  cost-effective  and  viable  for 
large-scale  production. 

Elemental  phosphorus  is  interesting  for  energy  conversion  and 
storage  purposes  since  it  possesses  distinctive  physical/chemical 
properties  and  electronic  structure  such  as  multiple  allotropies,  a 
number  of  oxidation  states,  and  so  forth.  Phosphorous  and  its 
compounds  have  been  utilized  to  develop  novel  energy  materials 
with  improved  performance.  Plenty  of  phosphorous  alloys,  such  as 
Ti-P  [6],  V-P  [7,8],  Mn-P  [9,10],  Fe-P  [11-18],  Co-P  [19,20],  Ni-P 
[21,22],  Cu-P  [23,24],  Zn-P  [25],  Ga-P  [26],  Mo-P  [5],  In-P  [27], 
Sn-P  [28,29],  have  been  reported  as  the  promising  anode  materials 
for  high  performance  batteries.  Among  these  alloys,  Fe-P  alloys  are 
especially  interesting  because  of  their  unique  properties  like  ther¬ 
modynamic  stability,  anti-corrosion  ability,  high  theoretical  ca¬ 
pacity,  environmental  compatibility,  and  low  cost.  According  to  the 
Faradic  Law,  the  theoretical  capacity  of  FeP4  can  reach  from  596  to 
1787  mAh  g_1  when  the  lithium-ion  transferable  number  ranges 
from  1  to  3.  However,  the  large  volume  change  during  lithium  ion 
insertion/extraction,  electrolyte  decomposition,  synthesis  diffi¬ 
culty,  and  poor  conductivity  of  FeP4  prevent  its  practical  application 
[14,17  .  It  has  been  reported  that  nanostructured  materials  with 
specific  morphology  like  porous,  layered  or  tube-like  architectures 
may  alleviate  the  volume  change  and  improve  the  electrode  reac¬ 
tion  kinetics  [5,18,24  .  However,  to  the  best  of  our  knowledge,  no 
report  has  been  made  yet  for  the  studies  of  cycling  performance, 
rate  capability  and  impedance  behavior  of  high  phosphorus- 
content  Fe-P  alloys  (such  as  FeP4). 

Herein  we  report  a  nanoscaled  and  high  phosphorus-content 
FePy  alloy  which  can  be  feasibly  synthesized  through  the  mecha¬ 
nochemical  method  by  using  common  raw  materials  of  iron  and 
red  phosphorus  powder.  The  electrochemical  performance  of  this 
material  has  been  investigated  by  a  number  of  techniques  as  well. 

2.  Experimental  section 

2.2.  Synthesis  of  materials 

The  sample  was  facilely  synthesized  by  the  mechanochemical 
method  in  an  inert  atmosphere.  Iron  (99.9%,  <10  pm,  spherical,  Alfa 
Aesar)  and  red  phosphorus  (98.5%,  100  mesh,  Alfa  Aesar)  powders 
were  mixed  together  in  a  molar  ratio  of  1 :4  (Fe  to  P)  and  transferred 
to  a  stainless  steel  can  in  a  glovebox  full  of  helium  for  ball  milling. 
The  weight  ratio  of  stainless  steel  balls  to  raw  materials  is  7.3.  The 
can  was  sealed  and  ball  milled  for  15  h  (every  1  h  milling  with  a 
30  min  break),  on  a  high  energy  SPEX  8000D  dual  mixer/mill  (SPEX 
SamplePrep  LLC,  USA).  After  being  cooled  to  ambient  temperature, 
the  sample  was  collected  for  the  following  characterizations. 

2.2.  Characterization 

The  phase  composition  of  the  as-made  material  was  determined 
by  powder  X-ray  diffraction  (XRD).  The  data  was  collected  from  15 
to  55°,  with  a  step  size  of  0.02°  while  spinning  the  sample  to 
minimize  preferred  orientation  on  a  Scintag  XDS2000  diffractom¬ 
eter  equipped  with  a  Ge(Li)  solid  state  detector  and  a  Cu  Ka  sealed 
tube  (A  =  1.54  A).  Synchrotron  powder  XRD  data  was  collected  on 
the  beamline  X7B  at  National  Synchrotron  Light  Source  with  a 
wavelength  of  0.3196  A,  and  quantitatively  analyzed  by  the  Rietveld 
refinement  approach  using  the  GSAS/EXPGUI  package  [30,31  .  The 


molar  ratio  of  Fe  to  P  was  determined  by  a  Varian  Vista-MPX  axial 
Inductively  Coupled  Plasma-Optical  Emission  Spectroscopy  (ICP- 
OES).  The  sample  was  dissolved  in  a  hot  aqua  regia  solution  and 
diluted  with  high  purity  water  for  elemental  analysis  [15  .  The 
particle  size  and  morphology  were  characterized  using  a  field 
emission  scanning  electron  microscope  (SEM,  ZeissSupra-55).  The 
elemental  composition  was  also  analyzed  on  an  energy-dispersive 
X-ray  spectrometry  (EDX)  by  dispersing  the  sample  powders  on  a 
silicon  wafer  using  ethanol  as  the  solvent. 

2.3.  Electrochemical  evaluation 

The  as-synthesized  powder  sample  for  the  electrochemical 
characterization  was  mixed  well  with  conductive  carbon  black  and 
polyvinylidene  fluoride  (PVDF)  binder  in  a  weight  ratio  of  80:10:10. 
After  adding  N-Methyl-2-pyrrolidone  (NMP)  solvent,  the  mixture 
was  mechanically  homogenized  in  a  stainless  steel  can  to  form 
viscous  slurry  to  coat  on  a  cleaned  copper  foil  using  a  doctor  blade, 
and  the  cast  thickness  was  accurately  controlled  to  be  20  pm.  After 
being  dried  at  80  °C  for  12  h,  the  copper  foils  loaded  with  active 
materials  were  cut  into  1.2  cm2  wafers  which  were  further  dried  at 
-100  °C  under  vacuum  for  -12  h  before  using  as  the  working 
electrodes.  Pure  lithium  foil  (thickness  -23  pm,  Aldrich)  was  used 
as  the  counter  and  reference  electrodes.  The  electrochemical  per¬ 
formance  was  evaluated  using  CR2325  coin  cells,  which  were 
assembled  in  a  helium-filled  glovebox  by  sandwiching  a  Celgard 
2400  microporous  separator  between  the  working  electrode  and 
metal  lithium.  The  electrolyte  was  1.0  M  LiPF6  dissolved  in  a 
mixture  of  Ethylene  Carbonate  (EC)  and  Dimethyl  Carbonate  (DMC) 
( 1 : 1  in  volume).  Electrochemical  performances  of  the  above  devices 
were  determined  between  0.1  V  and  2.0  V  vs.  Li+/Li  with  constant 
current  charge/discharge  cycles  at  298  K  under  various  current 
rates.  Cyclic  voltammetry  (CV)  tests  were  conducted  in  different 
voltage  ranges  to  obtain  the  proper  working  potential,  and  at 
different  scanning  rates  to  investigate  the  reaction  process.  Elec¬ 
trochemical  Impedance  Spectroscopy  (EIS)  measurements  were 
carried  out  at  different  cycles  with  a  sinusoidal  excitation  voltage  of 
10  mV,  and  the  impedance  curves  were  fitted  using  Zsimpwin  and 
Zview  software  [32,33].  All  electrochemical  tests  were  carried  out 
on  a  Biologic  VMP2  multichannel  potentiostat. 

3.  Results  and  discussion 

3.2.  XRD  analysis 

The  powder  XRD  patterns  of  the  starting  materials  and  as- 
synthesized  FePy  are  shown  in  Fig.  1,  exhibiting  iron  is  crystalline 
with  an  obvious  peak  at  -45°  while  red  phosphorus  is  in  amor¬ 
phous  form.  Black  crystalline  powder  of  FePy  is  formed  after  high- 
energy  ball  milling  for  15  h.  Phase  identification  reveals  that  this 
product  is  composed  of  FeP4,  FeP2  and  FeP.  A  three-phase  Rietveld 
refinement  of  the  synchrotron  XRD  data  is  shown  in  Fig.  2,  and  the 
refinement  results  are  summarized  in  Table  1.  The  refinement  re¬ 
sults  show  that  the  weight  ratio  of  FeP4,  FeP2  and  FeP  in  this  ma¬ 
terial  is  49.9%,  17.2%  and  32.9%,  respectively,  suggesting  the  as- 
made  material  is  a  composite  phase  comprising  different  Fe-P  al¬ 
loys.  According  to  the  ICP-OES  analysis  result  (Table  2),  the  molar 
ratio  of  P  to  Fe  is  3.8  without  dilution  and  3.7  with  dilution,  and  the 
slight  excess  of  Fe  can  be  ascribed  to  the  loss  of  Fe  from  the  stainless 
steel  balls  and  can  during  the  ball  milling  process.  It  should  be 
pointed  out  that  although  we  applied  the  stoichiometric  ratio  for 
the  starting  materials  (i.e.  the  molar  ratio  of  Fe  to  P  is  1 :4),  it  turns 
out  the  final  product  is  still  a  three-phase-mixed  composite.  It  is 
known  that  the  mechanochemical  method  utilizes  mechanical 
energy  to  make  a  chemical  reaction  take  place,  and  the  reaction 
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Fig.  1.  XRD  patterns  of  the  starting  materials  and  as-synthesized  FePy  (Miller  indices 
and  JCPDS  numbers  of  the  standard  patterns  are  provided). 


process  is  sometimes  complex.  Some  literature  have  mentioned 
that  the  mechanochemical  method  can  easily  create  defects  (such 
as  dislocations)  to  afford  non-stoichiometric  compounds,  and  the 
reaction  process/mechanism  can  be  material-dependent  [7,34-36]. 
Among  iron  phosphide  alloys,  to  the  best  of  our  knowledge,  FeP4 
single  phase  has  not  been  synthesized  through  the  mechano¬ 
chemical  method  yet,  as  other  impurity  phases  such  as  FeP  and/or 
FeP2  are  always  found  [14,37,38  .  This  might  be  due  to  the  partial 
sublimation  of  phosphorus  during  the  high  energy  ball-milling,  as 
well  as  the  special  thermodynamic  equilibrium  achieved  by  this 
specific  reaction. 

3.2.  Morphology  analysis 

According  to  the  SEM  images  in  Fig.  3,  the  FePy  material  consists 
of  particles  with  different  size.  A  large  number  of  small  particles 
(quasi-spherical  shape)  locate  on  the  surface  of  big  particles;  these 


Table  1 


Composition,  unit  cell  parameters  and  refinement  results  of  FePy. 


Chemical  formula 

FeP y 

Composition 

FeP4 

FeP2 

FeP 

Crystal  system 

Orthorhombic 

Orthorhombic 

Orthorhombic 

Space  group 

C222i 

Pnnm 

Pna2i 

a  (A) 

5.024(3) 

5.338(2) 

5.173(3) 

b  (A) 

9.870(2) 

5.820(4) 

5.797(4) 

c(A) 

5.668(3) 

2.959(4) 

3.150(2) 

V  (A3) 

281.081(2) 

91.932(2) 

94.438(3) 

0.0261 

Rwp 

0.0384 

Rf2 

0.0630 

big  particles  are  actually  aggregated  by  the  small  particles.  Due  to 
the  strong  collision  force  of  the  high  energy  milling,  the  size  of  the 
smallest  particle  is  even  below  20  nm,  which  might  be  beneficial 
for  mitigating  the  volume  change  during  lithium  insertion / 
extraction  and  thus  improve  the  electrochemical  performance  [5]. 
In  addition,  these  small  nanoparticles  will  supply  plenty  of  reac¬ 
tion  active  sites  for  lithium  ions  to  enhance  the  electrode  reaction 
activity.  According  to  the  EDX  results,  Fe  and  P  elements  are 
dominant  in  the  FePy  sample  (the  small  amount  of  Si  and  C  shown 
under  EDX  may  come  from  the  Si  base  and  carbon  coating  layer 
processed  for  SEM  test).  The  semi-quantitative  analysis  of  EDX 
shows  that  the  atomic  ratio  of  P:  Fe  is  -3.7,  matching  with  the  ICP- 
OES  result. 

3.3.  CV  analysis 

The  initial  open  circuit  voltage  (OCV)  of  the  fresh  coin  cell  (with 
FePy  as  the  active  electrode  material)  is  about  3.0  V.  CV  curves  at  a 
scanning  rate  of  0.01  mV  s-1  with  different  scanning  potential 
range  are  provided  in  Fig.  4(a)  and  (b).  The  reduction  peak  at  ~0.6  V 
is  coupled  with  the  oxidization  peak  at  -1.2  V,  forming  a  redox 
couple  for  lithium-ion  insertion/extraction  during  the  discharge/ 
charge  process.  Within  the  potential  range  between  0.01  and  2.5  V, 
the  redox  peak  intensity  decreases  gradually  during  cycling.  In 
addition,  no  scanning  area  appears  within  the  potential  range  of 
2.0-2.5  V,  indicating  that  the  capacity  contribution  from  this  range 
is  very  low.  Therefore,  we  changed  the  scanning  potential  range  to 
0.1 -2.0  V,  as  shown  in  Fig.  4(b).  From  this  figure  an  obscure 


Fig.  2.  Rietveld  refinements  of  synchrotron  X-ray  diffraction  data  of  FePy  to  identify  the  phase  composition.  Red  bar:  FeP2;  Green:  FeP;  Blue:  FeP4.  (For  interpretation  of  the  ref¬ 
erences  to  colour  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 
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Table  2 

ICP-OES  results  of  the  as-synthesized  FePy. 


Sample 

Fe  (ppm) 

P  (ppm) 

Fe/P  (in  mole) 

Blank 

0.01525 

0.00567 

Undiluted  sample 

6.54558 

13.70780 

1  :  3.8 

Diluted  sample  (1:10) 

0.66555 

1.34182 

1  :  3.7 

reduction  peak  can  be  observed  for  the  1st  cycle  while  the  oxi¬ 
dization  peak  position  is  slightly  lower  than  ~1.2  V  and  the  intensity 
is  not  very  high.  This  may  be  related  to  the  electrode  activation  as 
well  as  the  electrolyte  wetting  process.  The  activation  continues 
with  the  electrochemistry  progressing  while  the  electrolyte  grad¬ 
ually  soaking  into  the  electrode.  During  the  3rd  cycle  the  reduction 
peak  intensity  reaches  maximum  and  then  decreases,  whereas  the 
maximum  intensity  of  the  oxidization  peak  occurs  during  the  2nd 
cycle  (and  then  decreases).  The  decrease  of  redox  peak  intensity 
possibly  suggests  a  reaction  activity  decrease,  which  may  be  asso¬ 
ciated  with  the  capacity  fade  [33,39].  However,  the  decrease  of 
redox  peak  intensity  between  0.1  and  2.0  V  is  much  lower  than  that 
between  0.01  and  2.5  V,  therefore,  the  potential  range  of  0.1 -2.0  V 
was  chosen  for  the  following  tests. 

3.4.  Galvanostatic  discharge/charge  tests 

Because  Fe-P  alloyed  anode  has  a  large  volume  change  during 
the  discharge/charge  process  [17],  the  loading  density  of  the  elec¬ 
trode  material  will  affect  its  electrochemical  performance. 
Accordingly,  four  electrodes  with  different  loading  weights  were 
prepared  to  investigate  the  electrochemical  performance.  Consid¬ 
ering  the  activation  process  of  the  first  two  cycles,  the  3rd  galva¬ 
nostatic  discharge/charge  cycle  curves  (current  density: 
0.03  mA  cnrT2,  20  mA  g-1,  0.03C)  were  chosen  and  compared  in 
Fig.  5.  All  these  discharge/charge  curves  exhibit  similar  shape  and 
charge  potential  plateaus,  however,  the  capacity  and  discharge  in¬ 
termediate  potential  decrease  as  the  loading  density  of  active 


materials  increases.  This  indicates  a  bigger  polarization  associated 
with  a  lower  lithium-ion  transfer  rate  as  well  as  more  serious 
disturbance  from  the  volume  change  in  the  thicker  electrode. 
Therefore,  an  electrode  loading  density  of  -1.4  mg  cnrT2  was 
selected  to  obtain  the  optimized  results. 

The  discharge/charge  curves  of  FePy  from  the  first  seven  cycles 
are  shown  in  Fig.  6(a).  Except  for  the  1st  discharge  curve  possessing 
an  average  potential  plateau  of  -0.3  V,  all  other  curves  have  a  stable 
discharge  and  charge  plateau  at  -0.7  and  -1.1  V,  respectively. 
Moreover,  all  these  discharge/charge  curves  (except  for  the  1st 
discharge)  exhibit  a  similar  shape,  which  reveals  that  this  material 
has  a  similar  reversible  lithium-ion  intercalation/de-intercalation 
behavior  after  the  initial  lithium  intercalation.  During  the  1st 
discharge,  the  number  of  inserted  lithium  ions  is  -13.3,  which  is 
higher  than  12,  the  maximum  lithium  accommodating  number  of 
FeP4  according  to  the  proposed  reaction  mechanism  14].  However, 
the  number  of  removed  lithium  ions  decreases  to  -9.7  after  the  1st 
charge,  and  the  number  of  re-inserted  lithium  ions  is  similar  to  this 
value  after  the  2nd  discharge.  The  excess  discharge  capacity  during 
the  1st  cycle  may  be  ascribed  to  the  formation  of  solid  electrolyte 
interphase  (SEI)  film  that  usually  takes  place  below  1.0  V  vs.  Li+/Li 
for  anode  materials  17];  such  a  SEI  film  may  result  in  a  weak 
reduction  peak,  which  is  consistent  with  the  above  CV  results.  This 
irreversible  capacity  loss  results  in  a  low  Coulombic  efficiency 
(-73%)  for  the  1st  cycle.  The  subsequent  charge/discharge  curves 
change  little,  with  all  the  discharge  curves  having  two  obvious  knee 
points  at  -0.9  and  -0.5  V,  along  with  other  two  knee  points  at  -1.0 
and  -1.2  V  on  the  charge  curves.  Moreover,  it  can  be  observed  that 
all  the  discharge/charge  potential  plateaus  are  somewhat  inclined, 
(different  from  the  flat  potential  plateaus  of  LLfTisO^  [40]), 
implying  that  the  lithium-ion  insertion/extraction  process  of  this 
material  behaves  as  a  single  phase-type  reaction. 

Fig.  6(b)  shows  the  galvanostatic  discharge/charge  curves  of 
FePy  at  different  current  densities  from  the  3rd  cycle.  With  the 
current  density  increasing  from  0.03  to  0.3  mA  cm-2,  the  amount  of 
lithium  reacting  decreases  from  9.7  to  0.5,  and  it  returns  to  2.0 


Fig.  3.  SEM  images  and  EDX  analysis  of  the  as-synthesized  FePy. 
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Fig.  4.  CV  curves  of  FePy  at  a  scanning  rate  of  0.01  mV  s  1  with  different  scanning 
windows  (a)  0.01-2.5  V  and  (b)  0.1 -2.0  V. 


Fig.  5.  Comparison  of  electrode  loading  density  on  the  3rd  cycle  discharge/charge 
curves. 


Fig.  6.  Galvanostatic  discharge/charge  curves  of  the  as-synthesized  FePy  sample  at  (a) 
a  constant  current  density  of  0.03  mA  cm  2  from  the  first  seven  cycles;  and  (b)  current 
densities  of  0.03, 0.06, 0.15, 0.30,  and  then  back  to  0.03  mA  cm"2,  respectively,  from  the 
3rd  cycle. 

when  the  current  density  back  to  0.03  mA  cm-2.  This  shows  a  slow 
kinetic  behavior  of  this  material,  and  that  may  be  related  to  its  low 
conductivity  as  well  as  the  difficulty  of  building  the  new  structures. 

To  the  best  of  our  knowledge,  there  is  no  report  about  the 
cyclability  or  rate  performance  of  high  phosphorus-content  iron 
phosphide  anodes.  Therefore,  we  investigated  the  cycling  and  rate 
performance  of  FePy  at  different  current  densities  (Fig.  7).  As  we 
mentioned  above,  because  the  1st  discharge  capacity  of  this  ma¬ 
terial  always  contains  contributions  from  side  reactions  such  as  SEI 
formation,  we  chose  the  2nd  cycle  (stabilized)  capacity  as  the 
reference  point  for  the  cycling  retention  study.  As  shown  in  Fig  7(a), 
at  the  current  density  of  0.03  mA  cm-2  (0.03C),  the  reversible 
discharge  capacity  of  FePy  remains  -1089  mAh  g-1  after  10  cycles, 
approximately  73%  of  the  2nd  cycle  capacity  (55%  of  the  1st  cycle 
capacity).  When  the  current  density  increases  to  0.06  mA  cm-2, 
after  10  cycles  the  discharge  capacity  decreases  from  908  to 
581  mAh  g-1.  It  decreases  still  further  when  the  current  density  is 
increased  to  0.15  or  0.30  mA  cm-2,  but  the  cycling  behavior  be¬ 
comes  stable.  When  the  current  is  decreased  back  to  0.03  mA  cm-2 
the  capacity  increases  to  314  mAh  g-1.  The  coulombic  efficiency, 
except  -73%  for  the  1st  cycle,  remains  -97%  (at  the  current  density 
of  0.03  mA  cm-2)  or  even  higher  (at  other  current  densities  of 
0.06-0.30  mA  cm-2).  The  key  values  of  discharge  capacity  and 
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Current  density  (A  kg*1) 


Fig.  7.  (a)  Cycling  and  rate  performance  of  FePy  (at  difference  current  densities  of  0.03,  0.06,  0.15,  0.30,  and  then  back  to  0.03  mA  cm  2),  (b)  Ragone  plot,  and  (c)  improved  cycling 
performance  by  adding  FEC  as  the  electrolyte  additive. 


coulombic  efficiency  at  different  current  densities  are  summarized 
in  Table  3.  The  Ragone  plot  based  on  the  3rd  cycle  discharge  ca¬ 
pacities  (under  aforementioned  different  current  densities)  is 
shown  in  Fig.  7(b).  The  3rd  discharge  capacity  can  reach  up  to 
1436  Ah  kg-1  when  the  current  density  is  20  A  kg'1,  but  it  de¬ 
creases  to  62  Ah  kg-1  when  the  current  density  increases  to 
200  A  kg'1.  Therefore,  this  material  has  a  higher  energy  density 
than  graphene  (or  doped  graphene)  [41  ,  Li4Ti50i2  [42]  and  nano 


Table  3 


Summary  of  the  cycle  performance  of  FePy. 


Current  density  and 
the  cycle  number 

Discharge  capacity 
(mAh  g^1) 

Columbic 
efficiency  (%) 

0.03  mA  cm  2 

1st 

1984 

73 

2nd 

1486 

96 

10th 

1089 

97 

0.06  mA  cm  2 

2nd 

850 

99 

10th 

581 

98 

0.15  mA  cm  2 

2nd 

236 

99 

10th 

149 

98 

0.30  mA  cm  2 

2nd 

69 

98 

10th 

50 

99 

20th 

41 

99 

0.03  mA  cm  2 

2nd 

314 

97 

10th 

244 

99 

Sn02/V205  as  the  anodes  [43  ,  whereas  its  power  density  needs  to 
be  further  improved  for  future  applications. 

As  shown  in  Fig  7(a),  the  capacity  fading  of  FePy  needs  to  be 
reduced  in  order  to  get  a  better  cyclability.  It  is  known  that  addi¬ 
tives  such  as  fluoroethylene  carbonate  (FEC)  can  improve  the 
cycling  performance  of  LiM^CXi/graphite  cells  [44-46  .  Therefore, 
we  tested  the  cycling  performance  of  FePy  by  adding  10  wt.%  FEC  to 
the  pristine  electrolyte  of  1.0  M  LiPF6/EC:DMC  (1:1  in  volume),  and 
the  results  obtained  are  shown  in  Fig.  7(c).  After  20  cycles  at 
0.03  mA  cm-2,  the  capacity  of  FePy  is  -1000  mAh  g'1,  which  means 
the  capacity  retention  increases  -22%  by  adding  FEC  additive  to  the 
electrolyte,  much  better  than  the  one  without  FEC.  Also,  the  cycling 
stability  is  obviously  improved  even  though  the  1st  cycle  irre¬ 
versible  capacity  loss  increases  a  bit.  Such  an  improvement  of 
cyclability  may  be  attributed  to  adding  FEC  can  benefit  the  SEI 
formation  and  stabilization.  In  addition,  other  approaches,  such  as 
nanosizing  the  particle  along  with  morphology  control 
[5,7,18,24,47-49],  surface  modification  with  Sn/SnO  [16  ,  ZnO  [50] 
or  organic  carbon  21,51  ,  doping  [52  ,  alloying  53],  binder  modi¬ 
fication  [54],  could  also  improve  the  rate  performance  and  cycle 
stability. 

3.5.  EIS  analysis 

In  order  to  better  understand  the  electrode  reaction  and  ca¬ 
pacity  fade  mechanism  of  FePy,  EIS  measurements  were  carried 
out  on  the  electrodes  at  different  state  of  discharge:  one  fresh  cell 
with  others  discharging  to  1.0  V,  0.5  V,  and  0.1  V,  respectively. 
They  are  denoted  as  “1D1.0V”  (1st  cycle,  discharged  to  1.0  V), 
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Fig.  8.  (a)  EIS  curves  of  FePy  at  different  state  of  discharge,  (b)  EIS  fitting  results,  (c)  the  equivalent  circuit  model  without  SEI  component  and  (d)  the  equivalent  circuit  model  with 
SEI  component. 


“1D0.5V”  (1st  cycle,  discharged  to  0.5  V),  “1D0.1V”  (1st  cycle, 
discharged  to  0.1  V),  and  “cycled  D0.1V”  (66th  cycle,  discharged  to 
0.1  V),  respectively.  The  Nyquist  plots,  along  with  the  corre¬ 
sponding  equivalent  circuit  models  for  EIS  data  fitting,  are  pre¬ 
sented  in  Fig.  8.  Obviously,  their  shapes  are  greatly  different  from 
each  other.  For  the  fresh  cell  (with  OCV  of  3.0  V),  the  EIS  curve 
consists  of  a  depressed  semicircle  in  the  high-to-medium  fre¬ 
quency  region  and  an  oblique  line  with  a  slope  of  ~1  in  the  low 
frequency  region.  When  the  cell  is  discharged  to  1.0  V,  the  semi¬ 
circle  becomes  smaller,  indicating  the  interface  resistance  and 
charge  transfer  reaction  resistance  decrease  as  the  lithium  ions 
inserting  into  the  material.  When  discharging  to  0.5  V,  the  semi¬ 
circle  splits  into  two  depressed  ones,  which  become  more  obvious 
when  the  potential  reaches  0.1  V,  similar  to  that  of  the  discharged 
Si/C  composite  anode  [55  .  As  shown  in  Fig.  8(b),  after  65  cycles 
the  2nd  semicircle  becomes  much  weaker  despite  at  the  same 
discharge  potential  of  0.1  V,  indicative  of  a  different  electrode 
reaction  process  before  and  after  cycling. 

Considering  the  galvanostatic  discharge  curves,  the  discharge 
potential  plateaus  are  all  above  0.5  V  except  for  the  1st  cycle, 
suggesting  the  SEI  film  formed  below  1.0  V  [17  .  Accordingly,  two 
different  models  provided  in  Fig.  8(c)  and  (d)  were  adopted  to  fit 
the  experimental  curves  with  and  without  the  contribution  from 
the  SEI  resistance.  Taking  into  account  of  the  non-homogeneity 
such  as  porosity,  roughness  and  localized  distribution  in  the  sys¬ 
tem,  constant  phase  element  Q  is  used  to  represent  the  capacitance 
in  the  models,  and  a  depressed  semicircle  (Core-Element)  will 
appear  when  Q  is  placed  parallel  to  a  resistor.  The  intercept  of  the 
curve  on  the  Z'  axis  in  the  high-frequency  region  is  related  to  the 
ohmic  resistance  of  the  solution.  The  semicircles  in  the  high 


frequency  region  are  attributed  to  the  SEI  formation,  or  interfacial 
resistance,  or  these  two  combined  dependent  upon  the  discharge 
state  [55,56  .  The  semicircle  in  the  medium  frequency  region  is 
associated  with  the  charge  transfer  reaction  resistance,  and  the 
oblique  line  in  low  frequency  region  can  be  ascribed  to  Warburg 
impedance  which  is  considered  as  the  semi-infinite  chemical 
diffusion  of  lithium  ions  in  the  bulk  electrode. 

The  experimental  curves  can  be  fitted  well  by  using  the  models 
of  f^s(Qint^int)(Qct^ct)W  and  f^sQint(^int(QsEI^SEl))(Qct^ct)W  with 
ZsimpWin  and  Z-view  softwares,  and  the  fitting  parameters  are 
summarized  in  Table  4.  In  the  equivalent  circuit,  Rs,  Rmh  Qint,  RsEb 
Qsei.  Ret.  Oct.  and  W  refer  to  the  solution  Ohmic  resistance  (including 
the  contributions  from  electrolyte,  electrode,  current  collector, 
etc.),  the  interfacial  resistance  (between  active  material  and  other 
adjacent  substances  such  as  current  collector,  conductive  additive 
and  binder),  SEI  resistance,  the  constant  phase  element  of  the  SEI, 
the  constant  phase  element  of  the  interface,  the  charge  transfer 
reaction  resistance,  the  double-layer  capacitance  (around  active 
material  particles),  and  the  general  Warburg  impedance, 
respectively. 


Table  4 

Fitting  of  the  impedance  plots. 


Cell 

Rs  (Ohm) 

Rint  (Ohm) 

Rsei  (Ohm) 

Rct  (Ohm) 

W  (Ohm  sM/2) 

Fresh 

2.6 

1027.1 

— 

20,139 

561.7 

1D1.0  V 

2.1 

658.5 

— 

52.6 

117.9 

1D0.5  V 

1.9 

250.6 

68.6 

22.1 

21.6 

1D0.1  V 

1.3 

156.9 

110.7 

90.8 

10.3 

Cycled  D0.1  V 

5.1 

109.7 

8.7 

14.5 

22.3 
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As  for  the  EIS  models  of  intercalation  reaction  materials  like 
LiFeP04  [32]  or  LLfTisO^  [57],  the  Rct  and  W  share  one  Q.  because 
the  charge  transfer  reaction  takes  place  during  the  ion  diffusion 
process.  However,  the  W  is  here  connected  with  other  resistance 
elements  sequentially,  and  Rct  has  a  separate  constant  phase 
element  Q.  of  reaction  interface,  indicating  the  charge  transfer  re¬ 
action  and  the  ion  diffusion  take  place  separately.  This  is  because 
the  FePy  alloy  would  store  energy  through  a  conversion  reaction. 
For  the  conversion  reaction,  the  charge  transfer  reaction  usually 
takes  place  at  the  interface  between  active  material  and  electrolyte, 
independent  of  the  ion  diffusion  process. 

After  65  cycles,  the  semicircle  from  the  charge  transfer  reaction 
becomes  weak,  indicating  a  decrease  of  conversion  reaction  activ¬ 
ity.  Compared  to  the  “1D0.1V”  cell,  the  solution  ohmic  resistance  Rs 
and  Warburg  diffusion  impedance  W  of  the  “cycled  D0.1  V”  cell  have 
increased  a  lot,  suggesting  lithium-ion  diffusion  becomes  difficult 
after  cycling,  meanwhile,  the  interface  resistance  Rmh  SEI  film 
resistance  Rs ei  and  charge  transfer  reaction  resistance  Rct  have 
decreased,  indicating  a  reduced  reaction  depth  of  the  SEI  film  and 
electrode.  Therefore,  the  semicircle  decrease  after  65  cycles  might 
be  associated  with  the  SEI  film  breaking  down  during  cycling, 
which  has  been  demonstrated  to  be  significant  important  for  the 
cyclability  of  anode  materials  [58].  Therefore,  such  a  SEI  breakdown 
during  cycling  may  be  one  of  the  key  reasons  for  the  capacity  fading 
of  this  material. 


4.  Conclusions 

In  summary,  nanoscaled  FePy  (3  <  y  <  4)  material  has  been 
successfully  synthesized  by  a  facile  mechanochemical  method 
using  common  raw  materials  of  Fe  and  red  P  powders  in  an  inert 
atmosphere.  The  smallest  particle  size  is  below  20  nm,  and  small 
particles  aggregate  into  big  microsized  particles.  FePy  can  be 
cycled  reversibly,  which  may  be  due  to  its  nanoscaled  particles 
can  reduce  the  volume  change  during  the  discharge/charge 
process  and  stabilize  the  electrochemical  reaction.  Although 
there  is  a  capacity  loss  during  the  1st  cycle  associated  with  the 
electrode  activation  and  SEI  formation,  the  2nd  and  3rd 
discharge  capacities  can  reach  1486  and  1436  mAh  g-1,  respec¬ 
tively.  A  capacity  of  1089  mAh  g_1  (-73.3%  of  the  2nd  cycle  ca¬ 
pacity,  55%  of  the  1st  cycle  capacity)  can  be  delivered  after  10 
cycles  at  a  current  density  of  0.03  mA  cm-2,  with  the  corre¬ 
sponding  coulombic  efficiency  of  -97%.  By  adding  fluoroethylene 
carbonate  (FEC)  to  the  electrolyte,  the  cycling  performance  of 
this  material  can  be  improved.  Due  to  the  composition 
complexity  of  this  three-phase  composite,  as  well  as  the  amor¬ 
phous  feature  of  its  lithiated  products,  the  detailed  reaction 
mechanism  study  of  FePy  is  still  underway.  Actually  there  is  still  a 
debate  about  the  real  reaction  mechanism  of  high  phosphorus- 
content  FePy  (such  as  FeP4).  Does  it  completely  go  through  an 
intercalation  reaction  (FePy  +  xLi  -  LixFePy),  a  conversion  reac¬ 
tion  (FePy  +  3yLi  -  yLi3P  +  Fe),  or  these  two  combined?  A 
conclusive  answer  is  still  missing,  either  because  of  the  amor¬ 
phous  phases  formed  (which  XRD  does  not  see),  or  because  of  an 
unknown  lithiated  phase,  i.e.,  LixFePy,  which  was  proposed  based 
on  the  “similar”  compound  Li7MnP4,  while  the  structure  of  Lix_ 
FePy  is  still  unknown  [14  .  Therefore,  synchrotron-based  studies 
such  as  pair-distribution  function  (PDF)  analysis  (which  is  a 
particularly  powerful  tool  to  probe  the  amorphous  phase),  are 
needed  to  determine  the  detailed  reaction  mechanism  of  this 
material,  as  well  as  the  key  factors  that  govern  the  capacity 
fading.  This  will  help  to  improve  the  stability  and  rate  perfor¬ 
mance  of  high  phosphorus-content  Fe-P  alloys  as  the  new  anode 
materials  for  next  generation  of  lithium  ion  batteries. 
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